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ABSTRACT 

We calculate the two-loop quantum corrections, including the back- 
reaction of the Hawking radiation, to the one-loop effective metric in a 
unitary gauge quantization of the CGHS model of 2d dilaton gravity. The 
corresponding evaporating black hole solutions are analysed, and consis- 
tent semi-classical geometries appear in the weak-coupling region of the 
spacetime when the width of the matter pulse is larger then the short- 
distance cutoff. A consistent semi-classical geometry also appears in the 
limit of a shock-wave matter. The Hawking radiation flux receives non- 
thermal corrections such that it vanishes for late times and the total ra- 
diated mass is finite. There are no static remnants for matter pulses of 
finite width, although a BPP type static remnant appears in the shock- 
wave limit. Semi-classical geometries without curvature singularities can 
be obtained as well. Our results indicate that higher-order loop correc- 
tions can remove the singularities encountered in the one-loop solutions. 
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1. Introduction 



The work on two-dimensional (2d) dilaton gravity models has shown that they 
posses the features one is interested to understand in a realistic gravitational collapse 
(for a review see [|TJ). Especially interesting is the CGHS model ||, which is soluble 
classically and it is a renormalizable 2d field theory. This raises a hope that the 
quantum theory may be tractable, so that the properties of quantum black holes could 
be understood. As far as the quantum solvability of the CGHS model is concerned, 
one can find the physical Hilbert space of states, which is the matter Fock space 
i, 1, |, i 0, |. This exact result allows one to show that the quantum theory is 
unitary, because the dynamics is generated by a free-field matter Hamiltonian, which 
can be easily promoted into a Hermitian operator acting on the physical Hilbert space 
||. Furthermore, one can show that evaporating black hole geometries appear in the 
semi-classical limit of a unitary gauge quantization of the CGHS model || |8|, which 
shows that at least in 2d evaporating black holes can exist in a unitary quantum 
theory. 

Although the results of || || imply unitary evolution for the quantum CGHS 
black hole, one would like to see what is the end-state geometry (i.e. is it a remnant 
or the black hole completely evaporates with the information being returned through 
the Hawking radiation [1]). For this one needs the exact effective quantum metric, 
which can be obtained either from the exact effective action, or from the expectation 
value of the metric operator. So far only the one-loop perturbative approximations 
of the effective metric have been obtained || |9|, [K| |J, which are valid in the week- 



coupling region of the spacetime. Especially interesting is the BPP geometry [ITU 
which describes an evaporating black hole which ends up as a remnant. This same 
geometry arises at the one-loop level of the operator formalism || , which is consistent 
with the unitarity. One is then interested to see what is the effect of the higher loop 
corrections for the BPP solution, and therefore calculating the two-loop corrections 
is the simplest thing to do. In this paper we calculate the two loop corrections to 
the effective metric by using the operator formalism of M. The advantage of the 
operator formalism over the conventional effective action approach is that it gives 
the spacetime dependence of the effective metric automatically, while in the effective 
action approach one has to solve the effective equations of motion, which in the most 
cases cannot be done explicitly. Also at two loops there is a large number of diagrams 
one would have to evaluate in order to obtain the contribution to the effective action. 

In section 2 we review the operator formalism of ||. In section 3 we review the 
BPP solution in the context of the operator formalism, since it is a starting point for 
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our perturbative calculation. In section 4 we calculate the two-loop corrections. In 
section 5 we examine the two-loop semi-classical geometry. In section 6 we present 
our conclusions. 



2. The operator formalism 

We start from the classical CGHS action || 

S= f d 2 x^~g \e~* (R + (V0) 2 + 4A 2 ) - ±(V/) 2 1 , (2.1) 

where is a dilaton scalar field, / is a matter scalar field, g, R and V are determinant, 
curvature scalar and covariant derivative respectively, associated with a metric g^ v 
on a 2d manifold M. The topology of M is that of R x R. The equations of motion 
can be solved in the conformal gauge ds 2 = —e p dx + dx~ as 



e r = e 



\ 2 x + x--F + -F„ , / = / + (x + )+/_(x~) , (2.2) 



where 



F± = a± + b ± x ± + J X dy J V dzT±±{z) , (2.3) 
and T±± is the matter energy-momentum tensor 

T±± = \d ± fd ± f . (2.4) 

The residual conformal invariance can be fixed by a gauge choice p = <f>, and the 
independent integration constants are a + + a_ and h±. An equivalent form of the 
solution (2.2) which is suitable for our purposes, is given by 

F ± = a± + /3±x ± + f dy(x ± - y)T ±± (y) , (2.5) 

JA ± 

where 

M± M± 
a± = a± + J dyyT ±± (y) , b± = f3± - J dyT±±(y) . (2.6) 

It is clear from the solution (2.2) that the independent dynamical degree of free- 
dom is a free scalar field /. This conclusion also comes out from a reduced phase space 
analysis [BJ- Consequently the quantum theory is that of a free mass- less scalar quan- 
tum field /, propagating on a flat background ds 2 = —dx + dx~ with the dilaton and 
the conformal factor operators given by (2.2) || The matter energy-momentum 
tensor operator is defined as 



T±± = \ : d ± fd ± f 
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(2.7) 



The normal ordering in (2.7) is chosen to be with respect to creation and annihilation 
operators associated to coordinates x by 



/27T Jo ^/2uJk 
where uj k = \k\. 

The physical Hilbert space of the model is just a Fock space J-{dk) constructed 
from a\ acting on the vacuum |0). The model is unitary because the dynamics is 
generated by a free-field Hamiltonian 

H= / dkuj k a\a k + E , (2.9) 

J — oo 

which is a Hermitian operator acting on J 7 , where E is the vacuum energy. Con- 
sequently the states at t = \{x + + x~) = const, surfaces are related by a unitary 
transformation 

#(t 2 ) = e-^ (t2 -* l} ^(ti) . (2.10) 
One can also define the Heisenberg picture 

tf = e lHt ^(t) , A(t) = e lHt Ae~ lHt , (2.11) 

which relates the covariant quantization to the canonical quantization. For example, 



1 r°° dk 



f(t,x) = e tm f(x)e~ im = I -p= a k e i{kx - Wkt) + ale-^-^ , (2.12) 



where x = \{x + — x~). Similarly, the operator expressions (2.2) are the Heisenberg 
picture operators. 

Given a physical state \l/o; one can associate an effective metric to \l/(t) = e~ lHt ^Q 

as 

e ^ ff (t,x) = e Ht,x) |^ ; (2.13) 

where is the inverse operator of the Heisenberg operator (2.2). The geometry which 
is generated by e^ eff via ds 2 = —e^ eff dx + dx~ makes sense only in the regions of M 
where the quantum fluctuations are small. This will happen if 

yf\ (e 2 ^) - (e^) 2 | « (e*) . (2.14) 

The effective conformal factor e^ eff can be calculated perturbatively by using a series 
expansion |J 

oo 

(-\ 2 x + x~ -F)" 1 = e 4 ' (l-e* 5F)~ 1 = e' p0 Y / e n4>0 5F n , (2.15) 

n=0 
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where F is a c-number function, e ^° = — A 2 x + x — F and 5F = F — F . Then 

oo 

-A^+ar - F)- 1 ) = e*° ^ e n *° (6F n ) . (2.16) 

?1=0 

A convenient choice for Fq is 

F = (F + ) + (F_.) , (2.17) 
since then the lowest order metric is a one-loop semi-classical metric 

e -4> = -\i x + x - _ - (F_) . (2.18) 

is chosen such that it is as close as possible to a classical matter distribution 
fo(x + ) describing a left-moving pulse of matter. The corresponding classical metric 
is described by 

e-» = - X 2 x + A{x + ) - \ 2 x + x~ (2.19) 

A 

where 

M(x + ) = A f + dyyT° ++ (y) , A 2 A = f + dyT° ++ (y) (2.20) 

J— oo J— oo 

and T° + = ^d + f d + f . The geometry is that of a black hole of the mass 

M= lim M(x + ) , (2.21) 

£+— > + 00 

and the horizon is at 

x - = -A = - lim A(x + ) . (2.22) 

a;+^+oo 

In the limit of a shock-wave matter distribution, for which 

T° + = a5(x + - x+) , (2.23) 

we have 

M(x + ) = \ax+6{x + -x+) , A = . (2.24) 

A 

The asymptotically flat coordinates (r] + ,r]^) at the past null infinity are given by 

\ x + = e A " + , x~ = -Ae- A "" , (2.25) 

while the asymptotically flat coordinates {a + ,a ) at the future null infinity satisfy 

\ x + = e Xa+ , \{x~ + A) = -e- Xa ~ . (2.26) 

The corresponding Penrose diagram is given in Fig. 1. 
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Note that a change of coordinates —>■ ^ defines a new set of creation and 
annihilation operators through 

^=w,r^M m±+b ^ ■ <2 - 27) 

The old and the new creation and anhilation operators are related by a Bogoluibov 
transformation 

/oo 
dq{b q a qk + b\(3* qk ) , (2.28) 
-oo 

and the new vacuum is given by |0^) = S |0 X ). This allows us to take for \l/ a coherent 
state 







o: 



(2.29) 



where 1 0^) 



) is the vacuum for the coordinates (2.25), while 



.4 



dk[f (k)al k - f*(k)a- k ] 



(2.30) 



where fo(k) are the Fourier modes of fo(x + ). 



3. One-loop metric 

By using the operator formalism we can calculate perturbatively the effective 
metric (2.13) via the expansion (2.16). The lowest order semi-classical metric will be 
given by the expression (2.18) which requires a calculation of the expectation values 
of the T±± operators ||. One can show that 

(^^-^^(g) 2 , Wr-W = -j^ , (3-D 



where k = so that 

Z47T 1 



e -Po = e -<t>o 



2 

K i.^lx2„+ -I 1 F 7.4-/. .4- _+\(9fo\ 



1 \ 2 + - 

A xx 









(3.2) 

The expression (3.2) can also be obtained as a solution of the equations of motion of 
an effective one- loop action JlO | 

S eff = S -^J d 2 x^/=gRU ~ x R-kJ d 2 x^f^{R<\> - (V0) 2 ) , (3.3) 

where Sq is the CGHS action (2.1). 

Note that a natural choice for A± is A± = — oo. This choice makes the constant C 
infinite and b± — 0. We will ignore this infinity, since from the effective action point 
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of view C is a constant of integration whose value is determined from a requirement 
of having a consistent semi-classical geometry. By choosing C = ^/t[log(/t/4) — 1] 



one can obtain a consistent semi-classical geometry p0[ . In the case of the shock- 
wave matter this geometry is well defined in the x + > 0,x~ < quadrant. In the 
dilaton-vacuum sector (x + < Xq) the solution (3.2) becomes static 

e - P0 = e -0o = c _ \i x + x - _ * l g(-A 2 x+a;") , (3.4) 

and it is defined for o > a cr , where o = log(— \ 2 x + x~) is the static coordinate. 
At a = a cr there is a singularity, and this line is interpreted as a boundary of a 
strong coupling region. This is a common feature of the semi-classical metrics [Q 
, and a consistent geometry can be defined for a > a cr by imposing a reflecting 
boundary conditions at a = a cr . However, in the operator approach we do not 
impose reflecting boundary conditions. A consistent geometry is defined only in the 
regions of the spacetime where the metric fluctuations are small, and these coincide 
with the weak-coupling region e~^° > ||. Note that a curvature singularity occurs 
at a = a cr for C < ^fi;[log(fi:/4) — 1]. Hence this naked singularity will not appear for 
C > i«[bg(«/4) - 1]. 

For x + > Xq one obtains an evaporating black hole solution 



e 



-po — a -<Po — n I _ _ \2 



C + — - X 2 x + (x- + A) - -\og{-X 2 x + x-) . (3.5) 
A 4 



The corresponding Hawking radiation flux at the future null-infinity is determined in 
the operator formalism by evaluating 

(*o|T__(r)|*o> , (3.6) 

where T (£~) is normal ordered with respect to the asymptotically flat coordinates 

of the metric (3.5) at the future null-infinity || §]. The ^ coordinates turn out 
to be the same as the "out" coordinates (2.26) of the classical black hole solution, so 
that 

2vr(T__(r)) = ^[l-(l + AAe A -)- 2 ] . (3.7) 

The expression (3.7) corresponds to a thermal Hawking radiation, with Th = ^ 
O, |TT|. The Hawking radiation shrinks the apparent horizon of the solution (3.5), so 



that the apparent horizon line meets the curvature singularity in a finite proper time, 
at 



The curvature singularity then becomes naked for x + > xf. However, a static solution 
(3.4) of the form 

e -Po = e -0o = q_ x 2 x + {x- + A) - ^ log(-A 2 x + (x" + A)) (3.9) 

can be continuously matched to (3.5) along x~ = x^ if C = ^«[log(«;/4) — 1]. A small 
negative energy shock- wave emanates from that point, and for x~ > x^ the Hawking 
radiation stops, while the static geometry (3.9) has a null ADM mass. There is 
again a critical line a = a cr , corresponding to a singularity of the geometry (3.9). 
Note that the scalar curvature of (3.9) is bounded at x~ = x~ , and the singularity 
comes from the pathological behavior of e~^°, which becomes ill-defined for x~ > x~ . 
This singularity can be interpreted as the boundary of the region where higher order 
corrections become important. The spatial geometry of the remnant (3.9) is that of 
a semi- infinite throat, extending to the strong coupling region. The Penrose diagram 
of the one-loop geometry is given in Fig. 2. 

Note that in the operator approach the spacetime is always R 2 , with the geom- 
etry of Fig. 2 defined in the region where the metric fluctuations and higher order 
corrections are small. The one- loop approximation (2.18) is valid in the regions where 

e 2< t>°\(5F 2 )\ «1 . (3.10) 

The condition (3.10) will certainly break down for e~^° = 0, which also determines 
the position of the curvature singularities of the one-loop metric. Therefore one can 
see directly in the operator approach why the one-loop singularities represent a border 
of the strong coupling region. Also note that (3.10) is perturbatively equivalent to 
the condition for small quantum fluctuations (2.14), so that a well-defined geometry 
appears in the region e~^° > 0, which coincides with the weak-coupling region. This 
situation can be best described by the Kruskal diagram of Fig. 3. The region to the 
left of the curve ABCD in Fig. 3 corresponds to the strong-coupling region. 

The operator formalism can also resolve problems connected with a shift in the 
classical horizon due to quantum corrections. From (3.8) it appears as if the classical 
horizon x~ = —A has been shifted to a new position x~ = x^ < —A. One can 
calculate the corresponding Bogoliubov coefficients, and due to this shift, they will 
be described by incomplete Gamma functions in the late time approximation |12[ , in 



contrast to the ordinary Gamma functions when there is no shift [|Tl|]. As a result, the 
Bogoliubov coefficients will have a different asymptotic behavior for large frequencies 
from the standard case, and as a consequence the Hawking flux will diverge [[12] . 



However, at the one- loop order (T) is given by (3.7), which is finite. The way out of 



this paradox is provided by the fact that we are working in a quantum gravity theory, 
and therefore the space-time geometry fluctuates. Hence the position of the horizon 
fluctuates, and consistency requires that the effective horizon must be at x~ > —A. 



4. Two-loop corrections 

In order to calculate the two-loop corrections it will be useful to redefine F as || 

e~* = a + (3±x ± - X 2 x + x~ - [* dy(x - y) (0„| T ±± (y) \0 V ) 

Ja ± 

- \ f dy + (x + - y+) (T ±± (y) - (O v \T ±± (y) \0„)) 

JA ± 

= C - X 2 x + x- - ^\og\X 2 x + x-\- F + - F_ , (4.1) 
so that the new F± are given as 



dy(x ± - y)T ±± (y) , (4.2) 

where T = T — (0 V \ T \0 V ) and (0 V \ T \0 V ) = 0. Another convenient redefinition is to 
rescale f(x) to \/2irf(x). 

The left-moving sector gives corrections due to the matter quantum fluctuations 
H ||. We introduce the following operator ordering 

: T ++ (x 1 )T ++ (x 2 ) := T ++ (x 1 )T ++ (x 2 ) - (0 V \ T ++ {x 1 )T ++ {x 2 ) \%) , (4.3) 
so that 

(tt„| : f ++ ( Xl )f ++ (x 2 ) : |* ) = --L^^^^logl^xO - V (x 2 )\ , (4.4) 
where 

(%\f(xi)f(x 2 ) \0 V ) = — | log |77 + "(xi) -r] + (x 2 )\ - ^ log |?7 _ (^i) ~V~( X 2)\ ( 4 -5) 
has been used ||. 

The expression (4.4) is still divergent when x\ — > x 2 since 

d^dxjoglri+ix!) -T] + (x 2 )\ = — — + \D xl (rj) + 0{x 1 - x 2 ) , (4.6) 

(xi - x 2 ) 2 

where D x (j}) = |^ — § (f 1 ^) is a Schwarzian derivative. Expression (4.6) can be 
regularized by subtracting (x\ — x 2 )~ 2 from d\d 2 \og{rji — i] 2 ), which corresponds to 
using a new ordering || 

: T ++ (xi)T ++ (x 2 ) := T ++ (xi)T ++ (x 2 ) - (0 V \ T ++ ( Xl )T ++ (x 2 ) \0 V ) 

- (0 x \[A,T ++ (x 1 )][A,T ++ (x 2 )]\0 x ) . (4.7) 
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This ordering gives for 



2 ; ' 1 \df df 



1 r x _ 

~7T~5 II / dxi(x + - Xi ) (x 1 x 2 y l {\og(x 1 /x 2 )y 2 - 
Sir i=i"' A+ \ (^i — X V ) uxi ux 2 



(4.8) 



Now we will calculate (4.8) for the shock- wave matter distribution (2.23). In 
order to make the calculation simple and well defined, we take a matter pulse 

d+f (x + ) = ^ [6(x + - x+) - 6(x + - 4 - e)] (4.9) 

where e is the width of the pulse. When e — > (4.9) gives T ++ for the shock-wave. 
Let x + > Xq + e, then (^+) is given by 



/ / dx\dx 2 {x + — Xi)(x + — X2 



4e7r 2 A+ y x + ^^(log^i/a^)) 2 (xi-x 2 ) 2 / 

(4.10) 

By using the expansion (4.6), we obtain 



= -^ + -<) 2 (S^i±iM))! + (^ . (4.11) 

Since e is small, we have 

(^+)^-^(^ + -4) 2 ^ + -4) • (4.12) 

Note that in the limit e — > (4.12) will vanish, so that there is no two-loop correction 
in the right-moving sector for the shock-wave. In the case of arbitrary pulses 

^5F 2 + ) = C + (x + ) 2 + C' + x + + C'l (4.13) 

for x + > Xq + e. Note that 

C + = f ° ^ f° ^ dxxdxa (: T ++ ( Xl )T ++ (x 2 ) :) , (4.14) 

Jx Jx 

in accordance with (4.12). 

One can get a similar result by using a ^-function regularization. One starts from 
the formula 

/ dkke^'^ = —. - , 77! - 772 ^ , (4.15) 

Jo (771 - r] 2 y 

which can be used to rewrite (4.4) as 

1 r°° 

— J Q dkke tk ^-^d xlV d X2V d Xl fod X2 f , (4.16) 
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so that 



1 ^ fn~^~ poo 
5Fl) = — TT / dr}Ae x{r]+ - m) - 1) / dkke ik{m 

1 Sit 2 f i J-oo JO 



-»») 



i=l 

1 



8tt 2 

where 



dkk\F(k,r) + )\ 2 , (4.17) 



F(k,ri + ) = p + drie tkv (e x{ri+ - ri) - 1)^ . (4.18) 
j— 00 c/77 



By comparing (4.13) to (4.17), one obtains for the pulse (4.9) 

r°° dk 



a r 
4ir 2 e Jo 

where u = \og\(xQ + e) and v = logAxo . When e — > 0, we get 



2 

c<j~ 2 + z/~ 2 cos k(uo — v) 



(4.19) 



n f°° f 
C + = — dkk^-^ + 0(e 2 ) . (4.20) 
4?t 2 Jo u 2 {Xo ) 

The infinite integral / °° kdk can be replaced by a sum X^Li n , which can be regular- 
ized by using the ( function 

00 

C(s) = 5>- , Res>l , (4.21) 

n=l 

in the following way 

r-oo i-N _A/2 

/ kdk = lim / kdk = lim — 

JO N^ooJo N^oo 2 

/AT AT \ 

= & + 1 )-' JV ) = & S n -^ 1 

\n=l n=l / 

00 00 11 

= E«-iE 1 = C(-l)-k(0) = -^-7 • (4-22) 

n=l n=l iZ 4 

Therefore 

C + = -T^2 lo g _2 ( A 4) , (4-23) 

again a small negative constant, and lim e ^ C+ = 0. Note that a finite result in (4.23) 
has been obtained by analytical continuation of £(s) to s = — 1 and s = 0, which is 
the same as if an infinite number has been subtracted from (4.20). But this is what 
has been explicitly done in the regularization (4.8), so that it is no surprise that a 
similar result was obtained. 



For our purposes it will be instructive to regularize (4.19) when e is small but 
finite by a momentum cutoff. In that case 

r A dk . ek 
2AxJ 
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^ a dk . 2 ek ^ . . . t . 



When A is becoming large, we have 

where 7 is the Euler-Mascheroni constant. Therefore C + is negative if the width of 
the pulse is less then the short distance cutoff l c = Xxq/'jA and otherwise C + is 
positive. 

The right-moving sector determines the back-reaction of the Hawking radiation || 
[U . In this case we use a point-splitting method for regularizing the operator products. 
It amounts to calculating an appropriate limes (x2%-i — > xii) of the expression || 

2n 

(0,HI^r/(^i))|0,) • (4.26) 

i=l 

The expression (4.26) can be calculated by using Wick's theorem and by using the 
expression for a two-point function (4.5). A normal ordering can be defined by an 
appropriate subtraction of the products of the terms (x~ — xj)~ 2 and dkfk • • -difi 
from the expression (4.26) before taking the limes, such that one obtains a regular 
expression after taking the limes. A useful formula for doing this is (4.6). 

In the n = 2 case, one can define 

4vr 2 : T__(.x 1 )T__(s 2 ) := lim \im + — r^^-J^- 

13^11 14^12 ^OXi OX 3 OX2 OX4 [Xi — Xs) OX2 OX4 

| x 1 df_ df_ + 1 1 1 \ 

2 (x 2 — X4) 2 dx\ <9x 3 4 (xi — x 3 ) 2 (x 2 — X4) 2 ' 



where Xi denotes x i . This ordering gives 

1 , _ . _ ,, N 2 

32^2 

(T__(xx)> <T__(x 2 )> • (4.28) 



(0^| : T__(xi)T__(x 2 ) : |0„) = (d xl d X2 log 1 77 - rj (x 2 ) 



The expression (4.28) is still divergent when x\ — > x 2 , and it will be the source of 
divergence in 

(5F 2 ) = J J dx 1 dx 2 (x~ - xi)(x~ - x 2 ) [d xl d X2 log \r]~(xi) - r]'(x 2 )\) . 

(4.29) 

One way to regularize (4.28) is by changing the definition (4.27) || ||. This 
amounts to 

If _ ^ 2 

(0 V \ : T__(xi)T__(x 2 ) : |0^) = — — j d Xl d X2 \og\r] (x x ) -7] (x 2 ) \ - 



"' " 32tt 2 V 1 2 6 " " J ' V n (xi - x 2 ) 2 ) 
+ (r„(x!)> (T„(x 2 )> , (4.30) 
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which is finite for x\ — > x 2 due to (4.6). Consequently (5F_J is given by 

3^nJ A _ dx i( x ~ ~ x i) ^xi^ 2 logh~(^i) -*T(s2)l ~ ^ -x 2 ) 2 ) ' ^ 4 " 31 ^ 

which is finite for A_ finite. However, it is natural to take A_ = — oo. In that case 
the integral (4.31) is independent of x~ , which can be seen by making a substitution 

Xi = x~e u (4.32) 

so that 

. . 1 roo roo / _ (ti+tg) \ 2 

^F 2 ) = — J J dt l( it 2 e tl+t2 {\ - e u ){\ - e* 2 ) (e - (e* 1 - e t2 y 2 ) . 

(4.33) 

The integral (4.33) is a divergent constant, and the way it arises is similar to the 
one- loop constant C. Hence we will consider (4.33) as a finite constant, whose value 
is determined from some self-consistency requirement, and therefore 



SFi) = C_ . (4.34) 

Note that the result (4.34) can be also obtained with the ^-function regularization. 
Namely, by making the substitution (4.32) in the integral (4.29) we obtain 

I /* OO f OO 

( SF -) = ^J Q I ^(l-e^Xl-e-^-^)- 4 . (4.35) 

By repeated differentiation of (4.14) with respect to rj we obtain 

I r°° 

- -y A = l dkk 3 e ik ^ , (4.36) 



so that 



where 



For large / we have 



! roo 

C_=lim -/ dkkKMk)] 2 (4.37) 

«-oo 6 • 32tt 2 Jo 1 wi v ; 

JJk) = [ l dt(l - e _ *)e iW . (4.38) 
Jo 



so that 



^ ,. 1 r°° „,o/ k 2 + 2 2smkl 2 cos kl 

C_ = km / dkk 6 1 

Jo 



oo6-32tt 2 7o \k 2 (l + k 2 ) k(l + k 2 ) k 2 (l + k 2 )) 

1 / f°° ,, , k 2 + 2 2 r 00 , rr 2 sinx f°° , 2xcosa;\ 

= km - ^ - / dkk— — - / ax— — / ax— - 

i^oo 6 • 32tt 2 \Jo k 2 + 1 Z Jo Z 2 + x 2 Jo I 2 + x 2 J 

6-32Wo \ k 2 + lj ( J 
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The last integral in (4.40) can be rewritten as 

(rN roo Arp \ 

j£ + -exp(-e-^)) 

oo 

= lim [±iV(iV + 1) - h\ = £ n ~ h = C(-l) - h , (4-41) 

iV— >oo 

where 7 is the Euler-Mascheroni constant and we used 



/ e — ^ = -Ei(-e) = -7 - log e + 0(e) . (4.42) 

J 1 37 



Therefore we obtain 



a negative constant. 



5. Two- loop semi-classical geometry 

The analysis in the previous section implies that 

'SF 2 ) = C- + C + (x + - x+) 2 6(x + - x$) , (5.1) 



for narrow matter pulses. We have also shown that a regularization exists such that 
C_ is a negative constant, while C + depends on the shape of the matter pulse. Clearly 
the constants C± are regularization dependent. This regularization dependence is not 
a problem, since from the effective action approach point of view C± can be considered 
as integrals of motion (the two-loop equations of motion will be of higher order in 
spacetime derivatives, and therefore new integration constants will appear). Therefore 
the range of C± could be determined from some appropriate consistency condition. 
Also, in the following analysis we will show that the corresponding geometries are 
essentially determined by the sign of C + . Note that in the case of arbitrary pulses of 
compact support (SF 2 ) is quadratic in x + for x + outside of the support interval (see 
(4.13)), so that (5.1) is a good effective approximation in the general case. However, 
C+ can be then of both signs, positive or negative, depending on the matter pulse 
profile, which can be seen from (4.14) and (4.25). 

The two-loop metric can be written as 

ds 2 = -e^dx+dx- ; 



1 + e 2 *'(C_ + C + (x + - x+) 2 6(x + - x+))\ , (5.2) 
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where is the one-loop dilaton solution 



e -*o = c-a(x + -x+)d(x + -x+) - X 2 x + x~ - f log |A 2 x + x~| . (5.3) 

We introduce new constants a and (3 such that C_ = —a 2 and C + = ±f3 2 . As in the 
one-loop case, the relevant quadrant is x + > 0,x~ < 0. For x + < Xq the solution 
(5.2) is static, and it can be written as 



e^ 2 = e 9 



1 - e*a 2 



C + e a -%a , (5.4) 



where a = log(— \ 2 x + x~) is the static coordinate. The solution (5.4) describes a 
two-loop corrected dilaton vacuum. The corresponding scalar curvature will diverge 
on the curve 

The equation (5.5) has two solutions e~*° = ±a, which correspond to the one-loop 
singularity lines where C is replaced by C =F a. Since the shock- wave intersects first 
the C — a line (a > 0), this will be the critical line. The semi-classical geometry will 
be defined for a > a c , where C + e°" c — \a c = a. The curvature singularity will be 
absent for C > a — |(1 — log |). Therefore if we want to avoid a naked singularity at 
two loops, the BPP choice of the one-loop C has to be increased to 

C = a + |(log|-l) . (5.6) 

Note that in the quadrants x + x~ > naked singularities are present for any value of 
C. However, for the observer located at the right spatial infinity, these are located in 
the strong coupling region e^ < 0, where the two-loop approximation breaks down, 
and therefore can be ignored. Also note that for C_ = +a 2 no curvature singularities 
appear, and the only singularity comes from (f> becoming complex for e~^° < 0. 

For x + > Xq the solution (5.2) becomes 



e <f>2 = e 0O 



where 

= C + ^--\ 2 x + (x- +A)-%\og(-\ 2 x + x-) . (5.8) 
A 

It describes a two-loop corrected evaporating black hole geometry. The curvature 
singularity is given by the curve 
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The curve (5.9) can be parameterized as 



x + x = — e a 



x x - A 2 + C + [ } 

where C a = C + e CT — \a, and we have set A = 1. When C + = — (3 2 , the relevant 
branch of (5.10) is a curve which starts from (xq, —k/Axy) and it goes to x + = oo 
with an asymptote x~ = —A — (3. When C + = (3 2 , (5.9) is a closed curve, which 
starts from (xq , — k/Axg ), extends to x + = o;q + (aA//3)(A 2 + Z? 2 ) -2 ", where it turns 
back and ends up at x + = Xq line. 

The apparent horizon curve is given by the equation (9+02 = 0, which can be 
rewritten as 

s-lAl k 2C + (* + -* + )e-*7A» 

4A 2 x+ ~ e- 2 ^o - 3a 2 + 3C + (x+ - x+) 2 ' K ' 

When C + < 0, this curve starts from (x^, — A — k/Ax^ ), intersects x~ = — A — /?, and 

then asymptotically approaches this line as x + — > oo. This behavior implies that the 

apparent horizon must intersect the curvature singularity at some (xf,x~), so that 

for x + > xf there is a naked singularity. If C+ > 0, then (5.11) is a closed curve, 

which starts from (xq, — A — k/Ax^ ), extends to x + + 5/4/3, and it comes back 

to line. The naked singularity will appear in this case, but it will be located 

in the strong coupling region e~^° < 0, and therefore it can be ignored. Note that 

this does not happen in the C + < case, where the naked singularity lies outside the 

strong coupling region (see Figs 4 and 5). As we are going to see later, this behavior 

is correlated with the behavior of the Hawking flux. 

When (3 = 0, the equation (5.9) becomes a shock wave singularity equation (5.5), 
and the relevant root of that equation is e~^° — a = 0. This solution describes a 
singularity line of the one- loop metric with a smaller ADM mass C + — a. For 
a minimal allowed value of C (5.6), this becomes the one-loop metric with C = 
j (log \ — 1). For a very small or vanishing (3, the (3 terms can be neglected in 
(5.11) and one obtains the one-loop apparent horizon line equation. Therefore the 
intersection point of the apparent horizon line and the curvature singularity line for 
the shock- wave at two-loops is given by the one-loop expressions (3.8). The change 
in the intersection point (3.8) when (3 ^ can be evaluated perturbatively, and it is 
given by 

^=-£^=F*-+ <>(<*> 

= _ 2A 2 C + (xt - + 2 

na a 
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Therefore when C + < a naked singularity will appear in the region x~ > 
x~,x + > xf, unless we impose an appropriate boundary condition. In the shock- 
wave case (C+ = 0), one can impose a static solution (5.4) for x~ > x~ 



e H = e <Po 



1 - e*a 2 



, e -*> =C- X 2 x + (x~ + A) - £ \og(-X 2 x + (x- + A)) , 

(5.13) 

which can be continuously matched to (5.7) at x~ = x~ if 

<7 = a-£(l-log£) . (5.14) 

The metric (5.13) with the value of C given by (5.14) does not have a curvature sin- 
gularity at a — j log(— X 2 x + (x~ + A)) = a cr , and the naked singularity is removed. In 
this case a = o cr corresponds to a curve where e~^° is ill-defined (branch-point singu- 
larity), in a complete analogy with the one-loop case. However, when C + = —f3 2 ^ 0, 
there is no static two-loop dilaton vacuum solution which can be continuously matched 
to (5.7). The best one can do is to take 



e 4>2 = e <Po 



l-e 2 ^{a 2 + (3 2 {x + -x+) 2 )] , (5.15) 



for x~ > x~ and x + > xf, where is given by (5.13) and (5.14). Still, the 
naked singularity remains. Hence the two-loop corrections make the one-loop remnant 
unstable. Note that when C + = (3 2 , there is no curvature singularity for x + > 
Xq + (aA//3)(A 2 + /9 2 ) _ 2 j and therefore in this case there is no need for a sewing 
procedure. By examining the Hawking flux, we will see that this solution has a well 
defined flux of emitted particles at future null infinity. 

In order to calculate the Hawking flux we need the asymptotically flat coordinates 
(o- + , a~) at J^. These are given by 

\x + = e Xa+ , X(x~ + A) = -e~ xa , (5.16) 

where 

This can be written as 

a- = -j\og[-\(A + x-)] + ^±(A + x-)- 2 , (5.18) 

so that a non-zero correction appears at two loops. When C + < 0, the coordinate 
change (5.16) is well-defined for x~ < —A — ^, otherwise it is two- valued. When 
C + > 0, then this coordinate singularity is absent. A direct consequence of (5.18) is 
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that the Hawking flux at X + will not have the thermal form (3.7), which can be seen 
by evaluating (3.6) for £~ = o~ . This gives 

where the primes stand for derivatives with respect to o~ and x~ = —XAe~ x ' n . One 
then obtains 

9 lT \ , T , , 1 y V(~yA + ^A 2 ) - C + P 4 (y) + CjP 2 (y)] 
2n(T^)\ I+ = T(y) = - {y - A )*(C + + y*)* > ( 5 - 20 ) 

where y = x~ + A, 

P 2 (y) = -2y 2 + 3Ay - |A 2 , P 4 (y) = y 2 (-V + lOAy - 5A 2 ) , (5.21) 

and we have set A = 1. In Fig. 6 we give plots of T(y) versus y for C + < and 
C + > 0, respectively. The dotted line represents the classical (zero loops) Hawking 
flux 

1 AA 2 - yA 

When C + < 0, T(y) is close to T (y) for early times, but as y — > 0, i.e. for late 
times, T(y) goes to zero and then it diverges to —00 at y — —f3. Note that in this case 
there is a naked singularity in the weak coupling region. It has been also found at 
the one-loop order that naked singularities are accompanied by pathological Hawking 
fluxes H], and it has been argued there that such a catastrophic divergence of the 
flux is a generic feature whenever naked singularities occur. On the other hand, when 
C + > 0, the naked singularity is absent, and T(y) is finite and continuous in the 
semi-classical region. Moreover, it smoothly goes to zero for late times, so that the 
total radiated mass is finite 

M rad = / dy—- T{y) < 00 , (5.23) 
J-00 ay 

which are features not present at zero and one loop order. Actually, (5.23) is finite for 
the BPP solution, but this is done somewhat artificially, by taking (5.22) to be zero 
for y > + A. As a result of this discontinuous change in the flux, a shock-wave 
of negative energy emanates from (xf, x^~) (a thunder-pop), so that the total energy 
is conserved. The two-loop solution with C + > has a desired property that the 
Hawking radiation turns off itself for late times, but then becomes negative, and it 
goes to zero at y — (see Fig. 6). The appearance of this negative energy flux is 
a characteristic of a situation where all of the infalling matter gets out, since then 
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from the energy conservation it follows that M ra( i = 0, which can be only achieved if 
T(y) becomes negative for late times. This behavior has been observed in the case 
of the one-loop BPP solution with reflecting boundary conditions, when M < M cr 
so that a black hole is not formed, and hence all of the infalling matter gets out (it 
appears at x + = +00, the right future null-infinity) |L3] . In our case, the black hole 
forms and there are no reflecting boundary conditions, so that the infalling matter 
will reach the left future null- infinity (x~ = +00). This can be interpreted as matter 
getting out, provided there are no singularities at x~ = +00. However, the line 
x~ = +00 lays in the strong-coupling region, and without the knowledge of the full 
non-perturbative solution we can not say what exactly happens there. Still, all this 
indicates that the Hawking flux will deviate significantly from thermality for very late 
times. Also the behavior of T(y) is consistent with the creation of particle-antiparticle 
pairs, where particles reaching infinity give rise to positive T(y), while antiparticles 
carrying negative energy give rise to negative T(y) fT3|j . 

Note that (5.23) is only a function of A, C + = [3 2 and A, and it does not depend 
on the incoming mass M. Also M ra d(ft,A) diverges for small /3's, so that for a 
sufficiently small (3 one will have M ra d > M, and these solutions will violate the 
energy conservation. This can only mean that the higher order corrections have not 
been taken into account. On the other hand, if one wants to have two-loop solutions 
which are consistent with the energy conservation, f3 should be larger (which means 
wider matter pulses, see (4.25)) so that M ra d(/3, A) < M. Solutions with M ra d > 
are consistent with the appearance of a remnant, while the M ra d = solution is 
consistent with the complete evaporation of the black hole. When (3 = 0, then one 
has a BPP type solution, with a massless remnant. In that case 

M rad = M-k^>M , (5.24) 

and the energy conservation is insured by emission of a thunder-pop, of negative 
energy which appears because the remnant metric and the black hole metric 

are not sewn up smoothly. 

6. Conclusions 

The two-loop corrections to the effective metric take a simple form (5.1), which 
is valid for arbitrary matter pulses when x + is outside of the pulse. The constants C± 
are regularization dependent, and one can consider them as integrals of motion. The 
corresponding semi-classical geometries are essentially determined by the sign of C + . 
When C + > 0, consistent semi-classical geometries appear, where consistency means 
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absence of naked singularities in the regions of the spacetime where the semi-classical 
approximation is valid and the total radiated mass is less then the infalling mass. 
The solutions with C + = (shock-wave geometry) are special, in the sense that their 
geometry is essentially the same as the one-loop geometry, with the two-loop corrected 
static remnant appearing as the end-state. The Hawking flux is exactly the same as 
the one-loop flux, and a thunder-pop appears. However, when C + > 0, i.e. when 
the pulse has a reasonable width, the Hawking flux receives two-loop non-thermal 
corrections, and it becomes a continuous function of time. There is no remnant in 
the weak-coupling region, and the Hawking flux becomes negative for very late times 
(approaching zero) indicating particle-antiparticle pair creation. For sufficiently wide 
pulses one can have M ra d < M, so that the total energy is conserved, in agreement 
with the unitarity. Solutions with M rac i > are consistent with a remnant appearing 
in the strong-coupling region. However, in order to check this, one would need the 
full nonperturbative solution, especially because M ra d could then vanish, implying 
that all infalling matter gets out and the black hole completely evaporates. Note that 
appearance of remnants is somewhat unnatural in models where infalling matter does 
not couple to the gravitational field, and in order to insure the energy conservation 
one has to prevent the infalling matter to reach the future infinity. Apart from simply 
terminating the evolution by hand after some time done in the BPP case) , a 

more natural solution would be that a singularity in geometry appears at x~ = +oo. 
Also note that when C + > and C_ > no curvature singularities appear at all, 
although singularities in remain (0 becomes complex in the strong-coupling region). 
Whether this can be interpreted in favor of the no-remnant scenario it is difficult to 
say, since at this moment we do not know what happens in the strong coupling 
region e~^ a < 0. Still, it is very indicative that the two-loop corrections can cure the 
problems encountered in the lower loop approximations. This gives us a hope that 
the full non-perturbative solution will be well defined in the semi-classical regions of 
the spacetime (i.e. regions where the metric fluctuations are small), so that one can 
obtain a definite answer about the fate of the 2d black hole. 

Note that the back-reaction effect on the Hawking radiation flux is such that 
the Hawking radiation must deviate significantly from thermality in the late stages 
of evaporation. This behavior is expected from general arguments, and it is very 
encouraging that it is recovered in a concrete model. An analysis of the Bogolibov 



coefficients confirms this ||12|| , and it will be interesting to better understand the effect 
of the horizon fluctuations, which is described by the operator (2.22). 

When C- < 0, a naked singularity can appear in the dilaton vacuum sector if 
the one-loop constant C < a + |(log § + 1). However, this singularity does not affect 
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the Hawking flux since it is independent of C_ and C. Also the constant C does 
not depend on the incoming matter state, and hence it can be chosen freely, such 
that the naked singularity is absent. However, C + depends on the matter state (see 
(4.15)), and for narrow matter pulses C + is negative, resulting in the appearance of 
the naked singularity and the pathological behavior of the Hawking flux. The authors 



of |14] have conjectured that the back-reaction prevents naked singularities to form. 
The C + < solution seem to be a counterexample to this conjecture. However, 
one has neglected the higher order corrections, so it is still possible that in the full 
non-perturbative solution the naked singularities do not appear in the semi-classical 
regions. Also it is very indicative from (4.25) that C + is positive for reasonable matter 
pulses, i.e. pulses which are wider than the short-distance cutoff l c , which can play 
the role of a 2d Planck length. 
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Figure 1: Penrose diagram of the classical CGHS collapse geometry. 
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D 




Figure 4: Kruskal diagram of the two-loop semi-classical geometry when C+ < 0. 
The straight dashed lines passing through D and G are x~ = —A and x~ = —A — (3 
lines, respectively. The curve BD is the strong-coupling border e~^° = 0, and the 
apparent horizon curve FCG intersects the curvature singularity curve BCG in the 
weak-coupling semi-classical region e~^° > 0. 
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Figure 5: Kruskal diagram of the two-loop semi-classical geometry when C + > 0. 
The apparent horizon curve FC intersects the curvature singularity curve BC inside 
the strong-coupling region, bordered by the curve BD. 
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